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ABSTRACT 


Two prediction models for graphical prognosis are developed and 
evaluated. The primary feature of these models which represents a de- 
parture from previous graphical methods is the inclusion of the ter- 
rain-induced vertical velocities. These vertical motions may in turn 
produce vorticity changes, and thus height changes. 

The first model presented is comprised of two separate parts; the 
first part is a prediction of the 500-mb surface and the second part a 
prediction of the 1000-mb to 500-mb thickness. The two parts may be com- 
bined to produce a surface (1000-mb) prognosis. This model failed to pro- 
duce satisfactory results, and hence a second model was developed. 

The second model presented follows somewhat the development of the 
first, but avoids an assumption relating the surface and 500-mb winds 
that is believed to have contributed largely to the original failure. 
This model operates on the current 1000-mb and 1000-mb to 500-mb thick- 
ness charts; and, utilizing a fraction of the 500-mb height change as 
predicted by any other established method, produces a surface (1000-mb) 
prognosis. The model was tested in two ways; one with the terrain factor 
included in the prediction equation and one omitting the terrain factor, 
Both methods produced satisfactory prognoses. However, inclusion of the 
terrain factor consistently produced the better results, 

The writer wishes to express his appreciation for the assistance 
and guidance of Professor George J. Haltiner of the U. S. Naval Postgrad- 


uate School in the preparation of this paper. 
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PREFACE 


Graphical methods of numerical prognosis have recently received 
widespread attention since the development by Fjortoft 1] of a predic- 
tion model for the 500-mb surface. Fjortoft provided a graphical tech- 
nique for integration of the barotropic vorticity equation. Despite the 
relative simplicity of the method and the rather restrictive assumptions 
inherent in its development, tests by the Air Weather Service of the 
United States Air Force [2 | indicate accuracy comparable to that achieved 
by a trained forecaster using conventional techniques. 

Following the success of Fjortoft's model, several other more in- 
volved graphical models have recently been devised. Although none as yet 
have been subjected to the degree of testing that Fjortoft's model has, 
preliminary results are encouraging. 

The graphical methods devised thus far appear to have been least 
successful in mountainous terrain areas. The reason for this probably 
lies in the simplified boundary condition, ®)J=0 at p,, which has been 
assumed in the previous models. 

It is the purpose of this paper to attempt to incorvorate terrain 
effects into graphical techniques of prognosis. The development here is 


somewhat similar to the treatments of Estoque [3] and Reed (i ° 
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CHAPTER I 
FIRST PREDICTION MODEL 
The first model is comprised of two parts; the prognosis of the 
500-mb surface and a separate prognosis of the 1000-mb to 500-mb thick- 
ness. The two prognoses may be combined for a resultant surface (1000-mb) 


prognosis. The development here parallels that of Estoque (sik 


1. Development of the prediction equation for the 500-mb surface. 
The "vertical velocity" is assumed to vary thusly: 

tay, psd kg (Vs H)(E)+ ems ye)|i- ($2) | on 
The first term on the right represents the contribution due to the oro- 
graphic effect and is assumed to vary linearly with pressure (p/p,)» 
from a maximim at Po to zero at the top of the atmosphere. k is a scalar 
multiplier with possible values Pen pineanecncen zero and unity. It is 
assumed to equal unity for this development and hence will not appear 
in subsequent equations in this chapter. The second term on the right 
is the contribution due to the non-orographic effects, and is assumed 
to vary parabolically from zero at D=P, to a maximum of GQ), at P=Pce 
Inherent in this distribution is the assumption that Pc (to be taken as 
S00mbs) is a level of non-divergence where there are no orographic effects. 
The non-orographic contribution assumed here varies slightly from Estoque's 
assumption of a sinusoidal variation of vertical velocity [3]. 


Differentiating (1) with respect to Te yields 


Bw _ - 2% P- Pe. 
> p = (Ys: YH) -  (e, -p. Po-P | (2) 





Evaluating (2) at Pe gives 


iG 7 
ey: RIE (V, ar H | (3) 
2 PLR, a. 


The following form of the vorticity equation is assumed: 


exa ; ow 
—— +V-V fa 5p (hs) 


This form of the vorticity equation neglects the vertical advection of 
ae: : ouy _ 

vorticity, the vortex tube term, and ie de Since 5 = Ve » Eq. 

(3) shows that the only divergence considered at 500mbs is due to ter- 

rain-induced vertical velocities, 


Combining (3) and () yields 
oF ev. g=- BE TH 6 


In order to obtain conservation equations utilizing only 500-mb 
data, the surface wind must be eliminated from Eq. (5). This will be 
done by assuming that Vj=aV, where a is a scalar multiplier. Then 
(5) may be written 


3 —afg 
SF vy. V Yq = a VV ER ° (6) 





Taking mean values of f and /? for the area to be prognosticated, 


(6) may be written 


oY agt 2 an 
St +v-V( Sa at P, H ) = 0 ~ (7) 


With the assumption of geostrophic vorticity, this becomes 


Sf 4+Vv.0(2V2 +f +2422 -H)=0 5m 





3 2 £ 2, 
Z + VV) Lye f+ 284-H) =o ls 


as f and agf2H/p, are invariant with time. 
Introducing the finite difference formla for vorticity into Eq. 
(8) gives 
L+V. V)2@-2) +4 + oat =O (9) 
4a? Te 
Multiplying through by factor fd°/lmég, which is invariant with time, 


and expanding, produces (a) 
a 


tid 4 afd? f 2 
2(- -¥) iy eae ly er 2 : H ai 
mG Hm af 
(b) (c) (d) 
yd pam 2) + vito opty. te’ 208 y= 0 
om Hg ~ $a? eee Ning= “P. 
Considering terms of Eq. (10) separately: 
(a) This term may be written o/art(z - Z), as other terms are 
invariant with time. 
(b) The factor lm@g/fd® may be taken outside of V operator, as 
variation is small compared with the variation of (Z - Z). 
(c) Following Fjortoft Iz], this term may be written as 
V.7u(g), where , 
Alla Z. 
Oe d sin? cos¥ 
J(Y) = még 


O 
(d) Taking mean value of f for area to be prognosticated, this 


df 


term may be written 
am 2 ae 2 





Rewriting (10) with above changes gives 





—2L— 37 
ie 2d 
= (2-2) +V-O((Z-De J) + Se H| =o (12) 
With the notation 7 
FP q* 
Beret dH , 
Hm’ f 


Eq. (11) becomes 
2 (7 -2e V-Y Z-Z+J(y) +B] = 0 
Bt = 4 


As B and J(y) are invariant with time, (12) may be written 


a (Z-Z+J) +B) +V- VE-Z +d) + B= 0 (13) 


or, in Jacobian forn, 
o/s T g 5 eae 
2 (Z-z+d+8) +2 (Zz Z+J)+B)=O, cw 


which is the desired prediction equation. However, in using the above 
equation, the advection would be accomplished with the current 500-mb 
wind, which obviously may vary considerably with time, especially if a 
prognostic interval of 12 or 2) hours is employed. In order to arrive 
at an advective field which is more conservative with time, the Jacob- 
ian identity J(X,Y)=J(X+Y, Y) may be employed. Applying this to (1)) 


gives 


= (Z-Z + Jew + B) +2 JZ+Je TB, 7-Z +Jcq) +B) =O, 19) 


Thus the advection is accomplished with a much more conservative field, 
(Z+ J(y) + B), Z being the only quantity that varies with time, and this 
much more slowly than Z. 


As a preliminary test of results of the model developed herein, 


mn 





Fig. 1 was prepared. An initial west-east zonal flow of 50 knots was 
assumed at 5O00mbs. The surface wind was assumed to be equal in direc- 
tion and $ the magnitude of the 500-mb wind (i.e.; a= 3). Eq. (15) 
was used and a prognostic chart for a time of six hours was prepared, 
using height changes as approximated by the first term of Fjortoft's 
series pais The resultant prognostic chart shows a trough in the lee 


of the Rocky Mountains, as would be expected under the circumstances. 


2. Testing of the prediction equation for the 500-mb surface. 

The first step in the employment of Eq. (15) was the construction 
of the Z chart. This chart was prepared by conventional graphical methods, 
using a grid distance d of 10°km as suggested by Estoque (3). Next, the 
contours of J(~) corresponding to a grid distance of 10°km were added 
graphically to the Z chart. Assuming a value of $ for a, the quantity B 
was evaluated as 1.5 x 107¢, assuming standard values of the meteorolog- 
ical parameters, neglecting variations of m°, and taking d as 103km, as 
above. (Variations of m° are small, as a conformal projection with stand- 
ard latitudes of 30N and 60N was used, and prognostication was confined 
to middle latitudes). A terrain contour chart of North America was mlti- 
plied by the above constant and isolines of the quantity B were drawn. 
The B chart was then added to the Z + J(~) chart to produce the advective 
field (Z + J(y) + B) . The Z chart was then subtracted from the advective 
field to produce the (Z - Z+ J(w) + B) field, or the conserved quantity 
to be advected. 

The above two charts were superimposed in register, and isolines of 
(Z - Z+ 3(~) + B) advected geostrophically with the (Z + J(y) + B) chart 


for the chosen prognostic interval of 12 hours. 


~ 








The prognostic 12-hr height changes may have been recovered by use 

of Fjortoft's series oF 
Az =A(2-z +d +B) + 2(2-Z te) +B) 

and compared with observed values of height changes, but for this test a 
more direct method of evaluation was used. The quantity (Z - Z + J(w~) + B) 
was advected for the 12-hr period and compared with the observed values 
of the same quantity as computed from the observed chart at the end of the 
prognostic period. 

The results of four tests using the above techniques and assumptions 
were generally poor. Indications were that in some isolated regions the 
inclusion of the terrain factor gave better results than the use of the 


prediction equation 
- (Z-z +J¢)) ar : J (Z +Jt¢) Z-Z +Ji@)= O. (16) 


(Except for the different vertical velocity distribution assumed, this 
equation is essentially the same as Estoque's model ely. However, these 
regions of improved results were not constant in space nor were they pre- 
dictable, and the areas of poorer results generally equalled or exceeded 
the improved regions, resulting in generally poorer prognoses by inclusion 
of the terrain factor than by ignoring it. 

It was reasoned that the choice of a= $, i.e., that the surface 
wind is equal in direction and 4 the magnitude of the 500-mb wind (in 
connection with the terrain factor), may have been largely responsible 
for the failure of the method. Therefore, another technique was develop- 
ed using a equal to either 3, zero, or -4, as determined in the follow- 


ing manner. The current 500-mb and 850-mb charts were superimposed, and 


areas delineated where the 500-mb contours were approximately the same 





direction as, at right angles to, or opposite in direction to the 850-mb 
contours. The 850-mb chart was selected as the chart most nearly represent- 
ing the surface in the majority of the mountainous areas. Then, in pro- 
ducing the (Z + J(¥) + B) chart, the contours of B as drawn for the first 
technique were added, considered’ zero, or subtracted in the areas indicat- 
ed. The subsequent steps in this technique were as before. 

Results of Haine the new technique were slightly better than for the 
assumption that a is everywhere $. However, improvements realized were in- 
sufficient to justify use of the method. A prime difficulty in applica- 
tion was the determination of appropriate Beites of a, which appeared to 
vary rapidly over small areas, especially in situations where closed lows 
were present at the 500-mb or 850-mb levels. An extensive search was made 
to locate a synoptic situation where the differences between 500-mb and 
850-mb wind fields were clear-cut, in an attempt to evaluate the technique 
objectively. Such a situation was located, and results of the testing in- 
dicated that whereas the sign of a was correctly chosen by the above 
technique, the magnitude of a appeared to vary considerably from the as- 
sumed values of $, zero, and ~s. AS a consequence, the prognostication 
was not appreciably better than that using Eq. (16), and indications were 
that even if the sign of a were predictable, the magnitude would vary to 
such a degree that prognostications ignoring the terrain factor would be 


more reliable, 


3. Conclusions as to prognosis of the 500-mb surface. 

No improvement over Estoque's prediction model [3] is realized by 
inclusion of terrain effects on vertical motion by methods developed in 
this paper. It is believed that the major reason for the failure lies in 
the difficulty of expressing the surface wind in terms of the 500-mb wind 

8 





in the 500-mb conservation equation. 


4. Development of the Becton equation for the 1000- to 500-mb thick- 
nesSe 
A detailed development of the prediction equation will not be pre- 
sented here. The development follows closely the methods of Section 1 of 
this chapter, so a listing of the major assumptions and the final equa- 
tion will be shown, | 
Assumptions: 
(a) The vertical velocity distribution is as described in Eq. (1). 
(b) Temperature changes are adiabatic. The following form of the 


adiabatic thermal equation is employed: 


3t ca) te soe (17) 


(c) Vr Si and Vo VF are relatively small and are neglected. 
(d) Similar assumptions as those used in progressing from Eq. (10) 
to Eq. (11) in Section 1 of this chapter. 


With the notation 


2 Fe 
__ gat d ZH ere otto. 
G = lem and Ky * opt ygm 


the prediction equation evolves as 
- (1 = hea +J(¢) +G) +V+ V (F alae + JS(g) a G ) aa O, (18) 


or, in Jacobian form , 


= A oe ee ee 
S (h-Kaht Jen + @) +3 J(zh Kh +Jey +G)=0. 5) 


Due to the failure of the 500-mb model to produce an acceptable prog- 


nosis, no attempt was made to test Eq. (19). Rather, a different develop- 





ment was attempted in which the assumption V,= aV was not necessary, as 
it 1s believed that this assumption was the major cause of failure. The 


new development is presented in Chap. IT. 
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CHAPTER IL 
SECOND PREDICTION MODEL 
This model is designed to be applied in conjunction with a separate 
500-mb prognosis, such prognosis to be made by any method the user's ex- 
perience deems the most accurate. Fjortoft's method BE, is suggested. The 


development herein parallels that of Reed [i]. 


1. Development of the prediction equation. 
The vertical velocity distribution of Eq. (1) is also assumed for 
this model. Differentiating (1) with respect to pressure yields 
OW {? a om 
; Pp $ yw 2 J 
7 Sees (P- Ps) 
Ge” > 


and evaluation at 1000mbs gives 


2 ?P P £2 (7 re : 21) 


Assuming that temperature changes are adiabatic, the following form 


om) 


of the adiabatic thermal equation is employed: 


£. 32 )t ae ve se +o =O (22) 


+ 


By integration from p, to Pos (22) becomes 
fo 


oh ae 
—4V5:Vha or Cyd, (23) 
oe ae 
Fs 
The surface wind may be used for the advection in lieu of a suitable 


mean wind for the layer by considering the wind hodograph to be linear 


and the thermal wind relationship to hold, i.e. : 








Evaluating the right side of (23) by use of (1) leads to 


ia p ? 
of w =~-o KGB (Vs+>VH = | 
| Gade BAAN WH) 5°) dp 


SL CU {f- $08) Tap 


—e indicated integration, and noting that Po= ePc» gives 


xf ede co. VHF |, . 


Beets 
roan _ (PP pe 
— -3(p- -P,) i te ) 


¥; 


5 


which finally reduces to 


v[odp = 2g 2 (Vs-V V H)- tT ETR Wm. (2h) 
Substituting from Eq. (2h) into Eqe (23) produces 


mM 


on 


Evaluating Eq. ()) at Py yields 


OJo Z 
yt Vee +4) = > = Fe). , (26) 





MWe 





assuming that Vg= Vo. This assumption is used throughout this develop- 


ment. Combining (26) and (21) produces 
om /2 z 
+ Vs ‘V(Y, +§)= a eee (27) 
4c a Ps : 
Multiplying (25) by factor 3f/o- Pe gives 


2ongeye per 2 2 
a op TNS Phe gong (v-TH) | ; (28) 


which, upon substituing in (27), yields 


2s Vs -V (3+) = - 228° (ys -VH) 











— 37) 2h yy ‘Vh+2o0ka 
T I er ari THe (29) 


Collecting advective terms and time derivatives, and again noting that 


a eRe (29) may be written 
2(¥, tegeh) = ~ V5" Vi¥+26 +=2h+t+ usaie A] (30) 


Making the geostrophic approximation for ae leads to 
Him | 
de (Z,-Z,) Te Pe hn so 


= Ymg Sim IK9t 2 
Vg: of CAS oe St eH om 


which, after multiplying by factor fa°/lmeg, yields 





ue 








= afd | | 
20°25) + Imig B ; : 


= Ms -V/(Z,- 2.) + Ha ich Sts hadtgy + MES 
- aos ae (32) 
) 


where 


‘9? d* sinY cosy 
Jc) * af 


ra*g 
following Fjortoft (1]. ea is saeietaion 
po UKE > * ae 
af m*> 42 Lim go Fe 


(32) may be written 


S [Br Etch tty) + F/=V. VEe. +Chascg) + F| 


as J() and F are invariant with time. This is the prediction equation, 


which may be written in Jacobian form as 


> = iT 
Ble zorcreden +F| = P Jfo,Brz.+chelo Flo 


However, if (33) were to be used in prognosis, the advection would 
be accomplished with the current 1000-mb chart, which may vary consider- 
ably over a prognostic period of 12 or 2); hours. In order to arrive ata 
more conservative advective field, the Jacobian identity 

JOY) = d&X+yY,Y) 


may be employed. Applying this to (33) leads to 
Ss Fe 
= \Z,-Z, +Ch +Je@) +F | = 
itis B/E +Ch+Jeo4 F Z,~-Z, +Ch +Joat F| ; (3h) 


Uy, 





which accomplishes the advection with a more stable field than Zoe 
To recover height changes, the local change in the advected quantity 


over the prognostic period may be written 
a Zo -Zeot Ch + Ji +F]=-A 9 
or 


A\ Z,-2.+Ch|=-A ; (38) 


as F and J(y) are invariant with time, With the relationship h= Ze - 2, 


(35) may be written 


A\Z,-Z.+ Cz -C Zo) =-A (36) 
Solving for A 2, yields 
| a3 
Az, 2--\a raz. +c AZ 5| (37) 
Cx 
Utilizing Fjortoft's method (1] to evaluate AZ, gives 
Dee =| A a (38) 
= % [A + CAE. 
Substituting (38) into (37) yields 
AZ =o lA += ee AZ. (39) 
eo C+ C+! 


Testing indicated that As 2A, so (39) is simplified to 


| a 
Z¥ = A Ee + — ae _ (1,0) 
which further reduces to 


Ngee One) (4 Oss gan (141) 
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For simplicity of application, Eq. (1) was further simplified to 
ae F O35 Ae. (42) 
which was the equation utilized. However, testing indicated that hq. 


(1) may have given better results than the simplified form of Eq. (h2). 


2. Testing of the prediction equation. 

The prognostic period of 12 hours was chosen to test the prediction 
equation. Four successive 12-hr intervals were prognosticated starting 
with the 0030Z chart of 1h November, 1955. Since the purpose of this 
paper is to attempt improvement in graphical prognosis by the inclusion 
of terrain effects, the prediction equation was applied in two different 
forms, First, a prognosis was made applying Eq. (3h), and a second prog- 


nosis using 


= Z-Z+Ch+ Jey) sr 3 J Z+Ch tJoy Zo tCht ge 


which is simply (31) neglecting vertical motion induced by terrain. With 
minor differences, Eq. (43) is the prediction equation derived by Reed Th. 
In applying (3), the following steps were carried out: 

(a) The current Z, chart was constructed by conventional graphical 
methods, using a grid distance of 600Knm, 

(b) The current thickness chart was multiplied by the constant C 
and added to the Z, chart. In practice, using a value of 0.75 for C in 
this step was found to give better results than the derived value of 0.9, 
and consequently this value was used. 

(c) Considering k= 1, the coefficient of H in the quantity F was 
evaluated to be 7.35 x 107¢. A terrain contour chart was multiplied by 
this factor to produce the quantity F, and isolines of F were added to 


contours of J(y). The chart of F + J(y) (Fig. 2) was added to contours 
16 











produced in step (b) to produce the advective field (Zo+ Ch+ J(y~) + F). 

(d) The current Z, chart was graphically subtracted from the ad- 
vective field to produce the conserved quantity (Z, - Zed Ghee (Gtk). 

(e) Charts (c) and (d) were superimposed, and isolines of the con- 
served quantity advected geostrophically with the advective field for the 
prognostic period of 12 hours. A difficulty arose in the advection process 
along the west coast of the United States and Canada due to the extreme 
gradients of the advective field caused by the high coastal mountain 
ranges. It was reasoned that where the terrain gradient was very large 
the horizontal wind would not be fully translated into upslope motion, 
and therefore result in lesser vertical velocities. Thus the scalar multi- 
plier k in Eq. (32) may be appreciably less than unity in regions of such 
strong terrain gradients. The amount of reduction of k was determined by 
cases tested. This reduction factor may be taken into account by appro- 
priately reducing the gradients of F where the gradient is very large. 

The final form of F is given in Fig. 2. 

(f) The advected isolines of (Z, - 2,+ Ch+J(¥) + F) were sub- 
tracted graphically from the initial isolines of the same quantity to 
produce isolines of A. 

(g) Half of the 12-hr 500-mb height change was added graphically 
to A to produce the ON Tap chart. In practice, any suitable method of prog- 
iosticating the 500-mb chart could be employed, but for this test the 
actual 500-mb height change was used in order that the prediction equa- 
tions, Eqs. (3l) and (3), might be more objectively compared. 

In applying Eq. (3), the steps were identical to the above, with 
the exception that in step (c) only the contours of J() were added. 


The results of the above tests are displayed in Table 1, which com- 
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TABLE 1 


iL 2 3 by 


Correlation coefficient : A 0.832 0.8h7 Os 710 0.790 


pe ey es Oe ee ee SS ee ee ie eS ee a eee ee ee SE oe ae Se ee 
= a oe ow of ay Oo De 68 Oo 28 ow 8 28 oo ow © oe es So 68 oe SE ee 68 Ss om oe OR Oo 68 oF 68 ow om 88 8 ew ee oe oe a oe 
72 OS ow of am 28 om ow ow OS OD SE ee oF 68 © ow EF ow oO oO o® oo 28 oO 6 68 ow 68 oO ow OS 28 28 2 ow oe 6 oe OP oe es os = ee 
om aw Se ow om OD OP ap OO ow ow om oe OD ow 8 om 8 om ow = ow © om ee om oe em oe ew om OS ow oe om GE 88 28 68 ow oe om om oO oe om om 6 a © oO Oo o® 6! 8 o® o® oe ow 6 So 8 ow ow ee oe ee oe 
= 2 ow o@ o@ oe ow a SE 8 ow 6 6 OP on SS oe 288 6 ee ow = 28 Se ee 68 ee 6 6 68 68 = oe 6 88 28 om om ee © 8 28 ot 2 ee oe oe 
oD aw am OE Ow em ow oP Ow on CF a 28 on Op et om 68 et ew om om © OF ow 6 a Ow ow 6 6 28 oe Ow 8 om om om oe 28 oe oO oe ow o® oe 





A : Prognostication including terrain factor; computations for 
area west of 96W, i.e., Rocky Mountain area, 


B : Prognostication excluding terrain factors; computations for 
area west of 96W, : 


C : Prognostication including terrain factor; computations for 
entire area covered by grid. 


D : Prognostication including terrain factor; computations for 
area east of 96W, 


1 : 12-hr prognostic period commencing 0030Z 1) November, 1955. 
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pares height changes as prognosticated by the two methods to the observed 
height changes. The statistical parameters were computed using a 8 point 
square grid, which enclosed the area bounded by points 28N, 12hwWs; 55N, 
136W; 53N, 65Ws and 27N, 79W. Figs. 3, h, and 5 illustrate height changes 
for the period 1230Z 15 November to 00302 16 November, 1955. Fig. 3 shows 
observed height changes, Fig. the height changes as prognosticated by 
Eqe (3h), and Fig. 5 the height changes as prognosticated by Eq. (3). 


The grid encompasses an area Slightly larger than the above illustrations, 


3. Conclusions. 

The inclusion of the fgenan factor in this model results in materi- 
ally improved prognoses. 

However, the process of making a prognosis by this method is rather 
lengthy, and whereas the inclusion of the terrain factor into the pre- 
diction equation adds no more steps to the routine, (F may be combined 
with J() on a permanent chart and used for each prognosis), the advec- 
tive procedure is nuch more delicate. The advective field has areas of 
very strong gradient with the isolines of the conserved quantity closely 
packed and almost parallel to the advective field. Consequently consider- 
able care mist be taken to insure accurate advection. 

In constructing the contours of F, terrain irregularities were 
smoothed to some extent. Some smoothing is mandatory, as a minutely 
accurate F chart is not reasonable to use when considering large scale 
flow, Also, it is evident that the smoother the F chart, the easier the 
graohical constructions will be. Further testing might indicate an op- 
timum balance between detail of the terrain features and representation 
for large scale flow, 

As noted in Sect. 2, a value of O.75 was used as a multiplier of h 
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Fig. 5 





23 





in producing the advective field rather than the derived value of 0.9. 
More extensive testing would doubtless indicate the optimum value of 
this multiplier. Reed [li] found that a value of 0.5 gave good results 


in the three cases tested with his model. 


2h, 





Ce 
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